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ABSTRACT 

We calculate the vector current correlation function for light valence quarks in the deconfined 
phase of QCD. The calculations have been performed in quenched lattice QCD at T ~ lA5Tc 
for four values of the lattice cut-off on lattices up to size 128"^ x 48. This allows to perform a 
continuum extrapolation of the correlation function in the Euclidean time interval 0.2 < tT < 0.5, 
which extends to the largest temporal separations possible at finite temperature, to better than 
1% accuracy. In this interval, at the value of the temperature investigated, we find that the 
vector correlation function never deviates from the free correlator for massless quarks by more 
than 9%. We also determine the first two non- vanishing thermal moments of the vector meson 
spectral function. The second thermal moment deviates by less than 7% from the free value. With 
these constraints, we then proceed to extract information on the spectral representation of the 
vector correlator and discuss resulting consequences for the electrical conductivity and the thermal 
dilepton rate in the plasma phase. 
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I. INTRODUCTION 



Many of the properties of dcconfined, strongly interacting matter are reflected in the structure 
of correlation functions of the vector current. Its spectral representation directly relates to the 
invariant mass spectrum of dileptons and photons. Depending on the invariant mass regime one 
can get information on in-medium properties of heavy [l| and light quark bound states. For low 
invariant masses one enters the hydrodynamic regime where spectral functions reflect transport 
properties of the thermal medium [3|-l8|. In the case of the vector correlation function and its 
spectral representation this is characterized by a transport coefficient, the electrical conductivity. 
Heavy ion experiments start getting sensitive to the relevant regime of low mass dileptons/photons 
of a few hundred MeV, i.e. invariant masses that are of the order of the scale set by the QCD 
transition temperature Q. 

Spectral functions in the vector channel, or equivalently cross sections for thermal dilepton 



production have been extensively studied in perturbation theory |ld.lll|. Non-perturbative aspects 
have been included through the hard thermal loop resummation scheme [l^ . A straightforward 
application of the latter, however, breaks down in the low invariant mass regime. This finds its 
manifestation in a too rapid rise of the spectral density which, in turn, leads to an infrared divergent 
Euclidean correlator of the vector current [isf . This low invariant mass regime has been examined 
in much detail lisl and it has been demonstrated that the spectral function at low energies 
will increase linearly with energy giving rise to a finite, non vanishing electrical conductivity of the 
quark giuon plasma. 

Lattice calculations of Euclidean time correlation functions can give access to spectral functions 
[l^- Extracting these directly, however, is an ill-posed problem and progress in this direction has 
been made only recently through the application of statistical tools like the Maximum Entropy 
Method (MEM) [17|. Subsequently this method has been applied to also analyze correlation 
functions of the vector current 18l42d|. 

In this paper we analyze the behavior of the vector correlation function at high temperature, 
i.e. in the deconfined phase of QCD. Our first goal is to generate precise data for this correlator 
which allow for its extrapolation to the continuum limit in a time interval extending to the largest 
temporal separations possible at finite T. From the large r behavior of the correlator we then 
collect evidence for non-perturbative structure of the corresponding spectral function at low energy. 
This forms the basis for the extraction of information on the spectral representation of the vector 
correlator by starting with an ansatz for the spectral function as well as by exploiting the Maximum 
Entropy Method (MEM). 

The paper is organized as follows. In the next section we collect basic formulas for vector cor- 
relation functions and their spectral representation which are related to the thermal dilepton rate 
as well as the electrical conductivity. Section III introduces thermal moments of vector correlation 
functions. In Section IV we present our numerical results on the vector correlation function calcu- 
lated at T ~ 1.45Tc. In Section V we discuss an analysis of these results based on fits to the data 
as well as on an analysis using the MEM. In Section VI we discuss consequences for the calculation 
of dilepton and photon rates as well as the electrical conductivity. We finally give our conclusions 
and an outlook on future studies in Section VII. 
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II. THERMAL VECTOR CORRELATION FUNCTION, DILEPTON RATE AND 

ELECTRICAL CONDUCTIVITY 

A. Correlation functions 

Let us start by summarizing some basic relations between the Euclidean correlation function of 
the vector current, 

J^(t, x) = q{T, x)jf,q{T, x) , (III) 

at non-zero temperature and its spectral representation on the one hand and thermal dilcpton 
rates and the electrical conductivity on the other hand. 

We will analyze Euclidean time two-point functions of the current at fixed momentum, 



G^Ar,p} = J G^.(t,x) e'f^ , (II.2) 

with 

G^,(r,f) - (J^(r,f)4(0,0)) . (II.3) 

Here (...) denotes thermal expectation values. 

The vector current receives contributions from connected as well as disconnected diagrams. 
However, contributions from the latter are parametrically small. In high temperature perturbation 
theory they occur at 0{g^ In 1/g) and give rise to the difference between quark number and isospin 



susceptibilities [2l|. This difference has been found to be small at all temperatures in the high 



temperature phase of QCD 2^, l23[. Moreover, in calculations of the correlation function of the 
electromagnetic current, 

Jem ^ Q/qHt, x)jf,q^{T, x) , (II.4) 

/ 

the contribution of disconnected diagrams vanishes for three degenerate quark flavors as it is 
proportional to the square of the sum of elementary charges Qf. The connected part, on the 
other hand is proportional to the sum of the square of the elementary charges of the quark flavor 
/, C'em ~ X^/Q/- the following we will limit our discussion to the connected part of the 
vector current, as it is generally done in the analysis of thermal observables deduced from the 
electromagnetic current. 

The current-current correlation functions can be represented in terms of an integral over spec- 
tral functions, /9^^(ci;,p, T), which are odd functions of the 4-momentum, i.e. Pfj^^{uj,p,T) = 
— cj, —p, T). We denote by pu the sum over the three space-space components of the spectral 
function and also introduce the vector spectral function pv = poo + pa- With this we obtain the 
corresponding correlation functions. 

In the following we will limit ourselves to an analysis of Euclidean time correlation functions at 
vanishing three-momentum. Therefore we will drop the second argument of the correlation and 
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spectral functions and also suppress the argument for the explicit temperature dependence, e.g. 
Gh(tT) = Gh{t,0,T). Occasionally we also find it convenient to express dimensionful quantities 
in appropriate units of temperature. For this purpose we will introduce dimensionless variables, 
e.g. w = cj/T, t = tT and Gh = Gr/T^- 

B. Spectral functions 

The integral over the 0*'' component of the vector current is the net number of quarks (q-q) in a 
given flavor channel, uq^t) ~ J d^xJo{T, x). As the net number of quarks is conserved, it also does 
not depend on Euclidean time r, i.e. no{T) — no(0). Consequently the corresponding correlation 
function, Gqo{tT), is constant in Euclidean time and its spectral representation is simply given by 
a (5-function, 

P00(^^) = -2TTXqUjS{0j) , (II. 6) 

where Xg is the quark number susceptibilitjQ 

X, = -^j A^x (Jo(r,f)4(0,0)) . (II.7) 

The time-time component of the correlation function thus is proportional to the quark number 
susceptibility, Goq{tT) = — x^T. The vector correlation functions Gu and Gy therefore differ only 
by a constant, 

GuiTT)^X<iT + Gv{TT) . (II.8) 

The spatial components of the vector spectral function, pultu), increase quadratically for large 
values of w. In the free field limit one obtains for massless quarks, 

plZii^) = -2nT^Lo5{u:) , 

pfj°<=(a;) 2-kT'^uj5{uj) + 7^ tanh(a;/4T) . (II.9) 

Note that in this limit the contributions of (5-functions in poo and pa cancel in the vector spectral 
function pv{^) = Poo('^) + Pui'-^)- This is no longer the case in an interacting theory at finite 
temperature. 

At high temperature and for large energies corrections to the free field behavior can be calcu- 
lated perturbatively; the vector spectral function can be deduced from the calculation of one loop 



corrections to the leading order results for the thermal dilepton rate [2J] . This yields 



puioj) (1 + —) u}'^ ta.nh{uj/4T) , a;/T>l. (11.10) 

ZTT V TT / 

At energies oj<T perturbative calculations as well as the resummation of certain subsets of diagrams 
(hard thermal loop (HTL) resummation [l3|) become complicated as several scales of order g"T 
start becoming important. In fact, the straightforward HTL-resummation 
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is known to lead 



^ As discussed above we only consider the connected part of the vector correlation functions. Within this approxi- 
mation Xq is actually not the quark number susceptibility but the isospin susceptibility [21^ . 
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to a infrared divergent spectral function, phtl('^) which also leads to a divergent vector 

correlation function 



13| 



mperature 

mm, 



The limit cj — > is sensitive to transport properties in the thermal medium and the spectral 
functions need to be linear in uj in order to give rise to a non- vanishing, finite transport coefficient; 
in this limit the spatial components of the vector spectral function yield the electrical conductivity 

^ = ^ lim PiM . (11.11) 

T 6 u;^0 ujT ^ ' 

In the free field, infinite temperature limit also the spatial part of the spectral function contains a 
(5-function at the origin. Different from the time-time component, where the 5-function is protected 
by current conservation, this (5-function is smeared out at finite ternperature and the low energy 
part of Pa is expected to be well described by a Breit-Wigner peak 

p^^iu^) = x,CBW^^^^^^ , (11.12) 

which yields a{T) /Cem = 2xgCBn//(3r). In the infinite temperature limit the width of the Breit- 
Wigner peak vanishes; at the same time cbw ~^ Xg ~^ a-nd consequently the electrical 
conductivity is infinite in the non-interacting case. 

In the high temperature regime the time-time component of the spectral function receives per- 
turbative corrections, refiecting the perturbative corrections to the quark number susceptibility. 
To leading order this gives 

Poo = -2ttT^u:5{uj) (l - i^9\T)\ . (11.13) 



C. Thermal dilepton and photon rates 

The vector spectral function is directly related to the thermal production rate of dilepton pairs 
with squared invariant mass oj^ — p^, 

dA^i+i- ^ ^ aim Pv(.'^,P,T) , . 

dud^p "™6^3 (w'-F')(e"/^-l) ' 
where aem is the electromagnetic fine structure constant. 

The vector spectral function at light-like 4-momentum yields the photon emission rate of a 
thermal medium, 

dRj _ aem Pvjuj = \p\,T) 

dV ~ ^""^ e-/T-l ■ ^"-^^^ 
The emission rate of soft photons, thus can be related to the electrical conductivity, 

III. MOMENTS OF THE VECTOR SPECTRAL FUNCTION 



In addition to the vector correlation function itself we will calculate its curvature at the largest 
Euclidean time separation accessible at non-zero temperature, i.e. at tT = 1/2. The curvature is 
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the second thermal moment of the spectral functions at vanishing momentum, 

1 f°° duj fuj\n 



" n\ d(rr)" 



tT=1/2 



/ -f-)" '"}^) . , H^u,V, (III.l) 
io 2ti\t) sinh(w/2r) ' - > V ^ 



where n is chosen to be even as all odd moments vanish. These thermal moments give the Taylor 
expansion coefficients for the correlation function expanded around the mid-point of the Euclidean 
time interval, i.e. around tT ~ 1/2, 

G«(rT) = ^Gg"M--rr . (IIL2) 

It is obvious from Eq. IIII.ll that in a calculation of the n-th moment the spectral function is 
weighted by a term, which has a pronounced maximum at tanh(w/2r) — Lo/2nT. The n-th order 
moment thus is most sensitive to the spectral function at to /T ~ 2n and we may expect that higher 
order moments will come close to results obtained with perturbative spectral functions. 

In the infinite temperature, free field limit the integral in Eg. III. 51 can be evaluated analytically 
and one obtains for massless quarks [25I, 



free, n f , A + COfi^ {2tT tT) C0s(2tTtT) \ 

S^'^^rT) = TT (1 - 2tT) . ^ '- + 2 ^ 



Gy V'-; \-\- sin3(27rrT) ^ " sin2(27rTr) , 

Glr^irT) = T3 + Gi^^'irT) . (III.3) 

From this we easily obtain the moments of the free spectral functions. The first three non- vanishing 
moments are given by 



r<{0)Jree ^ 2 (q) J^ee ^ Q{2)Jree ^ 2871^^3 „(4),/ree ^ V2Ak^^^ 



We note that G^""* = g[7 ' for all 71 > not only at an infinite but at all values of the temperature as 
the correlators Gu and Gy only differ by a constant (Eq. lII.8( ). However, while in the free field limit 
Gh^tT) contains a constant contribution which drops out in the calculation of higher moments 
(n > 0), this is no longer the case at finite values of the temperature: all thermal moments, G^"^ 
with n > 0, will be sensitive to the smeared 5-function contributing to pii{uj) although, as indicated 
above, we expect this contribution to become more and more suppressed in higher order moments. 
In the following we will analyze the ratio of Gh{tT) and the free field correlator GjJ^'^^{tT), 



as well as the ratio of mid-point subtracted correlation functions 



(111.5) 



G(J-'=^(tT)-G^°^'^'^^'^ 

- Jfc (1 + - {I -rT)\ ....) . (III.6) 

Here we used the notation i?^'™"* = G^^^ /G^^\ Note that the curvature of these ratios at the 
mid-point determines the deviation of ratios of thermal moments from the corresponding free field 
values. While the ratios of correlation functions differ in the H = ii and H = V channels due to the 
additional constant contributing to Gv{tT), this constant drops out in the subtracted correlation 
function, i.e., Ay(rT) = A,,(tT) and G^^'^ = G|"^ for n > 0. 
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/? 


T/Tc 


csw 


K 


Zv 


16 


6.872 


1.46 


1.4125 


0.13495 


0.829 


24 (I) 
24 (II) 


7.192 


1.42 


1.3673 


0.13431 
0.13440 


0.842 


32 


7.457 


1.45 


1.3389 


0.13390 


0.851 


48 


7.793 


1.43 


1.3104 


0.13340 


0.861 



TABLE I. Simulation parameters for the generation of gauge field configurations on lattices of size 



p 




mAwi/T[fi = 1/a] 


^Ms/nf^ = '2GeV] 




= 2GeV] [MeV] 


6.872 


0.13495 


0.0182(4) 


0.0243(1) 




9.5 


7.192 


0.13440 


0.0182(4) 


0.0237(1) 




9.3 


7.192 


0.13431 


0.0816(2) 


0.1062(2) 




41.6 


7.457 


0.13390 


0.0790(3) 


0.0989(4) 




38.7 


7.793 


0.13340 


0.0964(2) 


0.1117(2) 




43.7 



TABLE II. AWI and MS' quark masses in units of the temperature. For orientation we show in tire last 
column estimates for the quark masses in units of MeV. These values have been obtained using for the 
phase transition temperature in the SU(3) gauge theory Tc — 270 MeV. 



IV. COMPUTATIONAL DETAILS AND NUMERICAL RESULTS 



All our numerical results are obtained from an analysis of quenched QCD gauge field configu- 
rations generated with the standard SU(3) single plaquette Wilson gauge action [2^. On these 
gauge field configurations correlation functions have been calculated using a clover action with 
non-pcrturbatively determined clover coefficients csw a-nd a hopping parameter k chosen close to 



its critical value 



27 



28| 



The gauge couplings used for calculations at four different values of the cut-off (oT = 1/-^t) 
have been selected such that the temperature stays approximately constant as the cut-off is varied; 
T ~ 1.45rc. For this we used an ansatz originally suggested in Ref. |29[. It is known to give the 
variation of the lattice cut-off as function of the gauge coupling to better than 1% in the interval 
[5.6, 6.5] js^l- We added to this analysis new results for the critical coupling (3c{Nr) and the square 

and extrapolated the fit results to the regime of couplings relevant 



31^ 



root of the string tension ^/a 
for our analysis, i.e. /3 G [6.8,7.8]. 

By tuning the Wilson hopping parameter k in the fermion action we have chosen quark masses 
that are approximately constant for our four values of the cut-off. We have estimated the quark 
masses using the axial Ward identity (AWI) to calculate the so called AWI current quark mass, 
ITT- AW 1 1 for the different values of the cut-of f 1321. Here we used a non-perturbatively improved 
axial- vector current with coefficient ca from |27l . |33| . 

To compare the quark masses at a common scale we first converted to renormalization group 

and then rescaled in the 



invariant quark masses 3j] using non-perturbative coefficients 
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MS" scheme to the common scale /x = 2 GeV using the four-loop perturbative running from [37| . In 
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/? 


# conf 


16 


32 


6.872 


251 




48 


6.872 


229 




64 


6.872 


191 




128 


6.872 


191 


24 (I) 


128 


7.192 


340 


24 (II) 


128 


7.192 


156 


32 


128 


7.457 


255 


48 


128 


7.793 


451 



TABLE III. Number of configurations analyzed on lattices of size A^^ x Nt. All configurations are separated 
by 500 Monte Carlo iterations. 

both cases errors are calculated from a jackknife analysis on tuawi and do not include systematic 
errors for the conversion to the MS* scheme. 

All basic parameters that enter our simulations are summarized in Table HI The values for 
the quark masses used in our calculations are quoted in Table HIl Our central results are based 
on calculations performed with renormalized quark masses nriy^/T ~ 0.1. Through calculations 
performed with a factor 4 smaller quark mass we checked on lattices with temporal extent Nr = 24 
that results on the large distance behavior of correlation functions agree within errors (see Fig.[T|). 
This smaller quark mass has also been used for the calculations on our coarsest lattice corresponding 
to a temporal lattice size — 16, where we analyzed the sensitivity of our results to finite volume 
effects. We thus expect that quark mass effects play no role in all results presented here. 

We have generated the gauge field configurations by using an over-relaxed heat bath algorithm. 
Configurations have been stored after every 500 sweeps. We have checked at several large tem- 
poral separations that the correlation functions calculated on these configurations are statistically 
independent. We further have measured the plaquette every sweep to obtain an estimate for 
the integrated autocorrelation time of order 1. Note that at temperatures above deconfinement 
topologically non-trivial configurations die out quickly. This is reflected in the observation that 
the inspection of the scalar and pseudo-scalar correlator revealed less than 3 % so-called excep- 
tional configuration which had an autocorrelation time of 3000 sweeps at most. The statistics 
accumulated on different size lattices is collected in Table IIIII 

For the computation of the quark propagators we used a plain conjugate gradient inverter. 
The algorithm has a stable convergence behavior, however, we found it appropriate to set the 
convergence criterion to 10"^"^ for the squared norm of the residue. This value was selected by 
monitoring the nearly exponential decay of spatial correlation functions which range over many 
orders of magnitude due to the large spatial extent of the lattices and screening masses of the order 
of 2tiT. 

For the calculation of the time-time (Goo) ^md space-space {Go) components of the vector 
correlation function we use the local vector current (r, x) introduced in Eq. III. II Unlike point- 



Note that zero modes do not couple into the vector and axial vector correlation functions. 
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split vector currents, this local lattice current is not conserved at non-zero lattice spacing. As a 
consequence the time-time component Goq{tT) may not strictly be r- independent at finite values 
of the cut-off [39I . Furthermore, the current needs to be renormalized multiplicativcly. Wherever 
required we therefore use the renormalized vector current, 

J^*(t, x) = i2KZv)qiT, xh^q{T, x) , (IV.l) 

with Zv denoting the non-perturbatively determined renormalization constant also given in 
Table H 

However, we note that many results presented in the following are given in terms of ratios of 
correlation functions as well as mid-point subtracted correlation functions divided by the quark 
number susceptibility, Xg = Xg/T'^- These ratios are independent of any multiplicative renormal- 
ization. 

It is obvious from Eq. IIII.5I that the correlators Gy (tT) and Ga (tT) are identical up to a 
constant, GooiTT) = Xg'^'^ ■ It thus suffices to analyze one of them. In the following we will mainly 
concentrate on a discussion of Gv{tT). This is because a subtle cancellation occurs between the 
constant contribution XqT^ to Gv{tT) and the Breit-Wigner term in the spatial part {Gu{tT)). 
In the free field limit this cancellation is exact. It is the incomplete cancellation of these two 
contributions to Gy {tT) which gives rise to the weak dependence of Gy / Gy'^'^ on Euclidean time 
and contains information on transport coefficients. Of course, the constant contribution XgT^ is 
trivial when it comes to an analysis of the spectral representation of vector correlation functions. 
For the determination of the spectral function puiuj) the correlator Gii{TT) thus is much more 
relevant. For this reason we will also show some results on the spatial correlator. 



A. Correlation functions of the vector current 



The vector current correlation function was calculated previously, using the same non-perturbatively 
improved clover action, on lattices of size 64'^ x 16 A MEM analysis was used to determine 
thermal dilepton rates which suggested a suppression of the spectral weight relative to the free 
spectral function at small energies. Subsequently vector correlation functions at high temperature 
were analyzed using staggered fermion formulations 19 , 2^ and a modified kernel was used in the 



MEM analysis [20| . This modification of the MEM algorithm increased the sensitivity of MEM to 
the low energy structure of spectral functions and resolved part of the problems observed in [18 1. 
We will show below that it is also crucial in a MEM analysis to choose default models which allow 
for a linear slope of the spectral functions at small energies. 

Calculations of vector spectral functions using staggered fermions are more involved, as two spec- 
tral functions, corresponding to different parity channels, need to be determined simultaneously. 
These studies were performed using lattices with temporal extent up to Nr = 14 Q and Nr = 24 



20|, respectively. They were performed with unrenormalized currents and primarily aimed at a 



determination of the (unrenormalized) electrical conductivity. They led, however, to quite different 
results, cr/T ~ 7Ce™ and cr/T ~ (0.4 ± 0.1)Cem 0- 

Here we improve over these studies with staggered fermions as well as the analysis of the thermal 
dilepton rate performed with improved Wilson fermions 18| by using the latter approach and 
doubling the spatial lattice extent to = 128 as well as increasing the temporal lattice extent by 
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1 35 .Gv(xT)/G*;^^(tT) 




FIG. 1. The vector correlation function, GvijT), calculated on lattices of size NlxNr &tT ~ 1.45Tc. The 
correlation functions are normalized to the free vector correlation function Gy'^'^{tT) in the continuum. 
Shown are data for tT > 1/Nt only. Label I and II refer to data sets generated with two different values 
of the quark mass (see Tabled. 



up to a factor three. We have calculated the vector correlation function on lattices of size x Nr, 
with 32 < iV^ < 128 and Nr = 16, 24, 32 and 48. For Nr = 16 we have calculated Gh{tT) on 
lattices with spatial extent N^, ~ 32, 64, 96 and 128. For Nr ~ 24 we checked that the quark 
masses used in our calculations are indeed small enough on the scale of the temperature to be 
ignored in the analysis of our correlation functions. On the largest spatial lattice, N^ = 128, we 
performed calculations for four different values of the lattice cut-off by choosing Nr = 16, 24, 32 
and 48 and at the same time changing the value of the gauge coupling j3 such that the temperature 
is kept constant, T ~ 1.45Tc. 

The large range of spatial lattice sizes used in this calculation, 2 < N^/Nr < 8, allows to quantify 
finite volume effects at fixed values of the lattice cut-off, aT = 1/Nr- The large temporal lattice 
allows to reduce the lattice spacing at T ~ 1.45rc to about 0.01 fm. Its variation by a factor three 
gives us good control over lattice cut-off effects. As will become clear in the following, finite volume 
effects are well under control and the variation of the lattice cut-off by a factor three now allows 
for the extrapolation to the continuum limit. This removes systematic errors that were present in 
earlier calculations performed with the same discretization scheme used here 18|. 

As the dependence of Gh{tT) on Euclidean time r is very similar to that of the free vector 
correlation function, we eliminate the dominant exponential variation of Gh{tT) with tT and 
show in Fig. [1] the ratio Gy [jT) / Gy'^'^ [tT) introduced in Eq. lIII.Sl for all parameter sets analyzed 
by us. We note that finite volume effects increase somewhat with tT. but remain small even 
for the largest Euclidean time separation tT = 1/2. Cut-off effects show the opposite trend; 
they increase with decreasing tT. In fact, they strongly influence the behavior of Gv{tT) at 
short distances for the first 6 to 8 Euclidean time units, 1/Nr < tT < 8/Nr- It also is obvious 
from Fig. [T] that deviations from the free field correlation function change qualitatively when the 
lattice spacing is reduced. For A^,- > 24 the ratio Gv{tT)/Gy'''^{tT) increases close to tT = 1/2 
indicating that the ratio of thermal moments, Ry'^\ is smaller than the corresponding free field 
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0.85 - 



32,x16 
48,x16 
6Cx16 
128 x1 6 



tT 
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1.2 
1.15 

1.1 
1.05 
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0.95 
0.9 

0.85 
0.8 



Goo(TT)/r 



128^x16 
128,x24 
128,x32 
128 x48 



0.2 



0.4 



0.6 



FIG. 2. The time-time component of the vector correlation function, Goo(Tr) /T^, calculated at T ~ 1.4521:; 
versus Euclidean time tT. The left hand part of the figure shows the volume dependence of Goo(tT)/T^ 
for N-r = 16 and 32 < iV„ < 128. The right hand figure shows the cut-off dependence of Goo(r, T)/T^ for 
Na ~ 128 and 16 < Nt < 48. In both figures the horizontal line shows a fit to the data obtained on the 
largest lattice, 128^ x 48. 





16 


24 


32 


48 


oo 




0.882(10) 


0.895(16) 


0.890(14) 


0.895(8) 


0.897(3) 




1.273(4) 


1.214(2) 


1.207(1) 


1.193(1) 


1.189(13) 


Gv(l/2)/(x,G(.-(l/2)) 
Gv(l/4)/(x.G{r-(l/4)) 


1.276(14) 
1.333(9) 


1.234(10) 
1.235(6) 


1.226(10) 
1.213(4) 


1.216(8) 
1.202(5) 


1.211(9) 
1.190(7) 


G„:(l/2)/(x,Gf;'-^^(l/2)) 
Gu{l/4)/{x,G{r{l/4)) 


1.184(7) 
1.302(7) 


1.156(6) 
1.213(6) 


1.152(5) 
1.193(3) 


1.145(6) 
1.183(5) 


1.142(9) 
1.172(7) 



TABLE IV. Quark number susceptibility (Xq/T^)) the curvature of the vector correlation function at the 
mid-point (Gy ^ ) and several values of the vector correlation functions expressed in units of the correspond- 
ing free field values and normalized with the quark number susceptibility. Results are from calculations on 
lattices of size 128"^ x A^'^. The last column gives the continuum extrapolated results. The quark number 
susceptibility has been renormalized using the renormalization constants listed in Table |I] In all other 
ratios renormalization constants drop out. 



value, Ry'°^ < Ry'/lee- ^hc spatial part alone, Gii{TT)/G{['"^ (tT), on the other hand drops when 
approaching tT = 1/2 (see Fig.©. This suggests > R^^jl^^. 

Also in calculations of the time-time component of the vector correlation function we find that 
finite volume, as well as cut-ofF effects, are small. In Fig. [5] we show results for —Goa{TT)/T^ for 
all lattice sizes analyzed by us. The left hand part of Fig.[2]shows results obtained on lattices with 
different spatial extent at fixed lattice cut-ofF. Except for the smallest aspect ratio, Na/Nr = 2, 
the results agree within statistical errors of about 1%. The right hand part shows results obtained 
on our largest spatial lattice, N„ = 128, at four different values of the lattice cut-off. This shows 
that also cut-off effects are small in the calculation of the time-time component of the vector 
correlation function, i.e. the quark number susceptibility. Apparently the fact that we use a 
non-conserved local current does not significantly alter the spectral properties of the time-time 
correlator. Except for short distances it is to a good degree r-independent. We summarize results 
for Xq/T"^ = -Goo{tT)/T^ calculated on the 128^ x Nr lattices in Table|lVl 
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We note that at T ~ lAbTc the quark number susceptibihty, Xq/T'^-, is about 10% smaller than 
the free field value Xg' ^^^^/T^ = 1 which is in accordance with calculations performed with staggered 
fermions. 

In Fig. [3] we show separately results for the space-space component of the vector current corre- 
lation function, Gii{TT)/T^ , and the combined vector correlator, Gv{tT)/T^ . In both cases we 
normalized the correlation function to the corresponding free correlator, GjJ'^''{tT)/T^. Of course, 
ratios ioT H — V and H = ii are related through Eq. III. 81 and provide identical information on the 
vector spectral function representing these correlators. 

It is apparent that the short distance part of the vector correlation functions is strongly influenced 
by finite cut-off effects. To some extent one can eliminate cut-off effects by calculating the free field 
correlation functions on lattices with finite temporal extent and considering the ratio Gy /Gy^'^'^"'^ 
rather than the ratio obtained by normalizing with the free vector correlator in the continuum 
limit. A comparison between the lower and upper panels of Fig. [3] indeed shows a reduction of 
cut-off effects, although the short distance part tT<0.2 still deviates substantially from the free 
field values. However, in the entire Euclidean time interval, where cut-off effects arc under control 
and a continuum extrapolation is possible for tT e [0.2 : 0.5], GviTT) as well as Gii{TT) stay 
close to the corresponding free field correlation functions. We also note that in the continuum 
limit the ratios shown in Fig. [3] will all approach the same value in the limit tT — >■ as all 
correlation functions diverge for tT -> and the fact that Gy and Gu differ by a constant becomes 
unimportant. Moreover, note that the data shown in this figure have been divided by the quark 
number susceptibility, x^/T^, to avoid the usage of any renormalization constants. Multiplying 
the data shown in Fig. [3] by the quark number susceptibility given in Table ITVl gives the ratio of 
the vector correlation functions and the free vector correlator. These values are about 10% smaller 
than the data shown in Fig. [3] We therefore expect that in the continuum limit deviations of the 
vector correlation function as well as its spatial component from the corresponding free correlator 
remain smaller than 9% for all tT > 0.2. 



B. Continuum extrapolation of the vector correlation function 

For three values of Euclidean time, tT = 1/4, 3/8, 1/2, we have numerical results from calcu- 
lations at all four values of the lattice cut-off. In these cases we performed directly a continuum 
extrapolation of the correlation functions using data from simulations on lattices with tempo- 
ral extent Nr = 24, 32 and 48. Specifically, we extrapolated the ratio of correlation functions 
G h{tT) / {xqGy'^'^ {tT)) by using a quadratic ansatz in aT = l/Nr- At other values of tT we 
use spline interpolations of the data sets on the Nr = 24 and 32 lattices to perform continuum 
extrapolations at Euclidean time distances available on the ~ 48 lattices, i.e. tT = fc/48 for 
9 < fc < 24. Continuum extrapolations at a few selected Euclidean time separations are shown in 
Fig. m We find that results obtained from the continuum extrapolation of correlators normalized 
with the free lattice and continuum correlation functions, respectively, agree within errors. The 
resulting continuum extrapolation of the vector correlation function obtained in this way is shown 
in Fig. [5j Data for tT = 1/4 and 1/2 are also summarized in Table ITVl 

From Fig. [5] we conclude that the largest deviation of Gy(rT) from the free vector correlation 
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^•55 -T2Gii{TT)/[x,G!r(TT)] 




FIG. 3. The space-space component of the vector correlation function, Gu{tT) /T^ , calculated at 
T ~ on lattices of size 128^ X A^^ versus Euclidean time tT (left) and the total vector correlator 

Gv{tT)/T^ (right). Both correlators have been normalized by the continuum version of the corresponding 
free vector correlation function defined in Eq. IIII.3I (upper part) and the discretized version of the free 
vector correlation function calculated on infinite spatial volumes with fixed Nt (lower part), respectively. 



function occurs at tT = 1/2. Taking into account the normalization with Xg obtain from 
Table HVl 



Gy(l/2) 
G(.'-'=^(l/2) 
Gv(l/4) 
1/4) 



G 



free. 



= 1.086 ±0.008 , 



(0.982 ± 0.005) 



Gvil/2) 



(IV.2) 



where the second relation has been obtained from a jackknife analysis of the ratio Gv'(l/4)/Gy (1/2). 

We note that the increase of Gy {tT) / Gy '^'^ {tT) with tT is significant. It becomes apparent 
only for sufficiently small lattice spacing, i.e. for large Nt, in particular in the normalization with 
the free continuum vector correlator. The rise with tT is a direct indicator that the vector spectral 
function in the low-w region is different from the free case, and is part of the motivation for the fit 
ansatze discussed in Section V. 
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FIG. 4. The ratio of the vector correlation function, Gv{t,T), normalized with the quark number sus- 
ceptibility and the free vector correlation function calculated in the continuum (full symbols) and on 
lattices with temporal extent Nt (open symbols). Shown are results for five values of Euclidean time, 
Tr = 0.1875, 0.25, 0.3125, 0.375 and 0.5 (bottom to top), on lattices with temporal extent A^r = 16, 24, 32 
and 48. For distances larger than tT = 0.1875 data have been shifted in steps of 0.1 for better visibility. 
For tT = 0.1875, and 0.3125 spline interpolations have been used on the Nt = 24 lattice to estimate results 
at these Euclidean time separations. Note that the far most right data set, corresponding to = 16, has 
not been included in the extrapolation. 
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FIG. 5. Continuum extrapolation of the vector correlation function. 

C. Curvature of the vector current correlation function 

Here we want to discuss the calculation of the second and fourth thermal moment of the vector 
spectral function defined in Eq. IIII.ll From Eq. IIII.5I we expect that the second thermal moment 
will be closer to the free field value than the correlation function at the mid-point itself. 

We analyze the ratio of the mid-point subtracted correlation function, Ay(rT), introduced in 
Ea. lIIL6l Results for Ay{tT) calculated at the four different values of the lattice-cut-off are shown 
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in Fig. [B] Again wc perform spline interpolations of results obtained on lattices with temporal 
extent Nr ~ 24 and 32. These interpolated data together with the results obtained on the Nt = 48 
lattices are then extrapolated to the continuum limit taking into account corrections of 0{{aT)'^). 
The extrapolated data at distances fc/48 are shown also in Fig. |6l These extrapolated data have 
been fitted to a quartic polynomial as indicated by the Taylor expansion given in Eq. IIIL6I From 
this fit we obtain 

(-.(2) 

— = 1.067 ±0.012 . (IV.3) 

^(2), free ^ ' 

Fits of the vector correlation function and the mid-point subtracted correlator are, of course, 
correlated. We made use of this and also calculated the ratio of the second thermal moment and 
the correlation function at the mid-point on jackknife blocks . We find 

V _ (n oao 4- n ni o\^v 



(0.982 ±0.012)- 



(0) ' (Q)Jree ' 

M2),free 

^ = (1.043 ± 0.010) . . (iv.4) 

^(0) ^ ' MO), free ^ ' 

Although the statistical significance of the deviation of Gy^ /Gy^ from the free case is marginal, 
the difference is consistent with the behavior of G(tT) / G^^'^^ {tT) close to tT = 1/2 shown in 
Fig. m At large enough tT the latter has a positive slope in tT , which approaches zero from 
above when tT reaches 1/2, and a small negative curvature that is in the limit of tT going to 1/2 
proportional to Ry'^'' — ^y/^ee' shown in Eq. IIII.5I 

We will show in the following that this puts stringent bounds on the magnitude of any contribu- 
tion to the vector correlation function that may arise from a peak in the vector spectral function 
at small energies. We also tried to look at the fourth thermal moment, which could be obtained 
from the curvature of the fits shown in Fig. |6l However, at present our numerical results do not 
allow to draw a firm conclusion about its value. 



V. ANALYSIS OF THE VECTOR CURRENT CORRELATION FUNCTION 



In the previous section we have presented results for the vector current correlation function in 
quenched QCD at T ~ 1.45Tc. For Euchdean times tT e [0.2,0.5] the results could be safely 
extrapolated to the continuum limit. As the correlation functions Gv{tT) and Gii{TT) are closely 
related to each other, we now concentrate on a discussion of the latter. Let us summarize our basic 
findings for the Gii(rT): 

(i) The correlation function at tT = 1/2 is about 2% larger than the corresponding free field 
value, G,,(l/2)/G^'-'=^(l/2) = 1.024(8) 

(ii) The deviation from the free field value increases with decreasing Euclidean time. At tT = 1/4 
the ratio is G,^(l/4)/Gf;"'(l/4) = 1.051(7) 

(iii) The second moment of the vector spectral function deviates from the free field value by about 
7%, G^.'VgI.'^'^'^'" = Gf/G^^^'f'^' = 1.067(12) 
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FIG. 6. The mid-point subtracted vector correlation function normalized to the corresponding difTerence 
for the free vector correlation function. Shown is the ratio as defined in Eq. IIII.6I but normalized by the 
quark number susceptibility. The continuum extrapolation is discussed in the text. 

This suggests that the spectral function of the free theory is a good starting point for the 
description of the vector correlator. We will use the knowledge about the asymptotic behavior 
of the vector spectral function at low and high energies, presented in Section II, to construct 
ansatze for the spectral function at high temperature. This is done in two steps. In the next 
subsection we discuss a straightforward implementation of the low and high energy constraints 
by using a combination of the free continuum spectral function and a Breit-Wigner contribution. 
Subsequently we will generalize this ansatz and discuss a class of spectral functions for which the 
continuum contribution is cut off at some low energy value and present also results from a MEM 
analysis. We find that the most straightforward approach already gives an excellent description of 
all our data. 



A. Breit-Wigner + continuum ansatz 

We start with the analysis of our results for the vector correlation function by using the ansatz, 

^00(1^) = ~2TrxqUjS{uj) , (V.l) 

pu{io) = 2x,CBW ^^ +^r/2)2 ^^^^^ tanh(^/4r) . (V.2) 

This ansatz depends on four temperature dependent parameters; the quark number susceptibility 
Xq[T), the strength {cbw{T)) and width (r(T)) of the Breit-Wigner peak and the parameter k{T) 
that parametrizes deviations from a free spectral function at large energies. At high temperature 
and for large energies, uj/T ^ 1, we expect to find k{T) ~ ag/Tr. Note that k{T) will also depend 
on w and actually will vanish for w — > 00 at fixed T. Wc will treat here k(T) as a constant and will 
not take into account any running of . As will become clear in the following, already this minimal 
ansatz provides a good description of current numerical results for the vector correlation functions 
Gu{tT) and Goo{tT). In fact, it is quite straightforward to use the latter for a determination of 
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The parameters of the Breit-Wigiier term can be extracted from fits to the vector correlation 
function. As pointed out in Eq. 111.81 the correlation functions Gy and Ga agree up to an additive 
constant, the quark number susceptibility. It thus suffices to analyze one of them. 

With the ansatz for the spectral function given in Eg. IV. 21 and F = F/T we obtain 

Gu{tT) = (1 + fc(r)) G{;'%tT) + CBwXgFBwirT, F) , (V.3) 

with 

T f°° ,^ w cosh(w(TT- 1/2)) , 

Fbw[tT,T) = — Au^ \ \ ,^,J " . V.4 

^ ' 2nJo (F/2)2+S2 sinh(c^/2) ^ ^ 

With this normalization one has limp_^g Fbw{tT, F) = 1 and limp^^ Ebw^jT, F) = 0. 

Note that with this ansatz the dependence on the continuum contribution and thus the depen- 
dence on k{T) can easily be eliminated by taking appropriately weighted differences of Gh{tT) 
at two values of Euclidean time. Similarly we obtain for the mid-point subtracted correlation 
functions, introduced in Eq. IIII.61 

A ^ 1 _u j_ ~ fBM/(rr,f)-FBiy(l/2,f) 

AylrJ ) = 1 + kil ) + CBWKa — ^7 ^7 . (V.5) 

^ ' ^ ' G(r^"(Tr)-G(r^"(i/2) ^ ' 

At tT =1/2 this fit ansatz yields results for the zeroth and second moment of the spectral function 
as introduced in Eq. IIII.5I and Eq. Illl.Si 

G,,(l/2) = 2(1 + k{T)) + CBwXqF^Bw(^) , 

Ay(l/2) = 1 + fc(T) + CBWX.j^j^, , (V.6) 

Gy ' 

with 



(2n)! 27r ((F/2)2 + sinh((^/2) 



In the limit F the thermal moments F^j^^ vanish for all n > and the fit ansatz for Av(rr) 
becomes a constant which relates to the deviations of the vector correlation function from the free 
field correlator at short distances. For all F > 0, however, the right hand side of Eq. IV. 51 is a 
monotonically increasing function of tT. We note that even for our largest lattice this feature is 
not yet obvious from the data. However, it becomes apparent in the continuum extrapolation of 
the mid-point subtracted correlation function. Having established the monotonic rise of A\/(rT) 
is important for motivating our fit ansatz. 

Fitting the continuum extrapolated correlation function GuItT) in the interval [0.2 : 0.5] to- 

(2) 

gether with the second thermal moment Gy to constrain the fit, we obtain 

k = 0.0465(30) , f = 2.235(75) , 2cBwXg/f = 1.098(27) . (V.8) 



This three parameter fit has a /d.o.f- = 0.06 for 12 degrees of freedom. The small /d.o.f. 
clearly reflects that even after the continuum extrapolation data at different distances are strongly 
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FIG. 7. Data for the continuum extrapolation of T'^Gv(tT) / {xqGy'^'^ {tT)) and the fit result for fixed 
Cflw/r and k[T) (left). The three curves show the result from a fit in the interval tT € [0.2 : 0.5] (central) 
and results obtained by varying F within its error band. In the right hand figure we show the spectral 
function obtained from the fit and compare with the free spectral function. 



correlated. Nonetheless, the fit provides an excellent description of the data. To illustrate the 
sensitivity of our fit to the low energy Breit-Wigner contribution and its dependence on Euclidean 
time, we show the fit to the data for Gy(rT) normalized to the free vector correlation function 
and the quark number susceptibility in Fig. [T] The error band shown in this figure corresponds to 
the width of the Breit-Wigner peak. The spectral function obtained from this fit is shown in the 
right hand part of the figure. Here also the error band arising from a variation of the width T is 
shown. 

It is clear from Fig. [71 that the vector correlation function is sensitive to the low energy, Breit- 
Wigner contribution only for distances rT>0.25. Taking into account also the value of the second 
thermal moment, the fits to the large distance regime return fit parameters which are well con- 
strained. As a consequence we obtain a significant result for the electrical conductivity, which is 
directly proportional to the fit parameter csw/F, 

T^— iTo — r~ ^ ' ^^-^^ 

which (accidentally) is close to the result found in j2oi | using staggered fermions with unrenormal- 
ized currents. It is more than an order of magnitude larger than the electrical conductivity in a 
pion gas |4C | . 



It should be obvious that this determination of the electrical conductivity is sensitive to the 
ansatz made for the spectral function in our analysis of the correlation functions. With this simple 
ansatz we obtain good fits of the vector correlation function with a very small chi-square per degree 
of freedom. However other ansatze may provide an equally good description of the current set of 
data. We will explore this in the next subsection by generalizing the current ansatz. 

We also note that the value determined for the correction to the free field behavior at large 
energies k ~ 0.05 at T ~ 1.45Tc is quite reasonable. Using the relation to the perturbative result, 
k = as/TT yields for the temperature dependent running coupling g^{T) = Anas ~ 2 which is in 
good agreement with other determinations of temperature dependent running couplings at high 
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energies or short distances 



B. Breit-Wigner plus truncated continuum ansatz 

In the ansatz used in the previous subsection the continuum part, i.e. a contribution proportional 
to the free spectral function, contributes to puiuj) for all uj. Also in the highly non-perturbative low 
momentum region there thus is a contribution proportional to ufi jT arising from the continuum 
part. In order to analyze its influence on the low energy structure of the spectral function we 
smoothly truncate the continuum contribution at some energy wq- We thus replace the spectral 
part of our fit ansatz by 

= IxqCBw +^r/2)2 + i + tanh(w/4T)e(cjo, A„) 

e(wo, A^) = (l + e(-o--')/-^-) . (V.IO) 

In the limit A^j — > the function 0(^0, A^^) becomes a 0-function with discontinuity at wq- We 
also used other cut-off functions for the continuum contribution, which lead to similar conclusions. 
We prefer the above version as it insures that the continuum contribution vanishes exponentially 
at w = 0. This allows, for instance, to replace the simple free field ansatz, used for the continuum 
in Eg. lV.lOl by the HTL or further improved spectral functions, which have power-law divergences 
at small w. 

As discussed in the previous section we perform three parameter fits with cbw^ T and fc as 
free parameters for several values of and A^j. We find that fits become worse with increasing 
Wo and/or increasing A^^. In both cases eventually too much of the continuum part at high 
energies gets suppressed. For small values of and A^^ the Breit-Wigner term compensates for 
the continuum contribution that has been cut off by increasing the low energy contribution, i.e. 
the intercept at a; = (electrical conductivity) rises. Results from fits which all lead to /d.o.f. 
smaller than unity are shown in Fig. [51 As a-nd A^^ increase the ^ j d.o.j . of the fits shown in 
this figure rises from its minimal value of about 0.06, obtained for loq/T ~ A^/T = 0, to unity. 
All fit parameters corresponding to the curves shown in Fig. [S] are summarized in Table |V] 

In particular the second moment of the correlation function normalized by the correlation func- 
tion at the mid-point, i.e. the ratio R^'^^ introduced in Eg. IIII.51 reacts quite sensitive to the 
truncation of the continuum part of the spectral function. This is shown in Fig. [51 where we com- 
pare the ratio R\'^'^'^ extracted from our continuum extrapolated data (error band) with fit results 
obtained for different values of ujq. The dependence on A^, on the other hand, is less pronounced. 

We conclude from this analysis that for a;/T>(2 — 4) acceptable spectral functions should not 
deviate significantly from the perturbative or free field like behavior. Moreover, the structure of 
the spectral function in the energy range a;/T<2 is sensitive to the form of the fit ansatz. Within 
the class of functions analyzed here the fits do, however, favor a small value for the cut-off ojq and 
a small value for the intercept of /9ii(w)/w at w = 0. 

These findings are consistent with the default model dependence of a MEM analysis which we 
will briefly discuss in the next subsection. 
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FIG. 8. Spectral functions obtained from fits to the vector correlation function using the ansatz given in 
Eq. IV. 101 For comparison we also show only the continuum part of the spectral function. The left hand 
figure shows results for different values of the cut-off (cjo) and fixed width (A^^). The right hand figure 
shows results for fixed uq/T = 1.5 and several values of A„. The curve labeled 'cont' is the continuum 
contribution to the fit described in Eq. IV. 31 
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TABLE V. Parameters for the fits shown in Fig. [8] left (upper part) and right (lower part). The last column 
gives the /dof of these fits. 



C. Analysis of vector correlation functions using the Maximum Entropy Method 

So far we did not make use of the Maximum Entropy Method that has been used in most other 
lattice studies of the spectral function in the vector channel. Here we want to discuss to what 
extent the analysis presented in the previous subsections can be reproduced in a MEM analysis, 
or whether a MEM analysis may improve over the result obtained with an ansatz for the spectral 
function. We performed a MEM analysis of the renormalized vector correlation function on our 



finest lattice using the modified kernel which has been introduced in 20j and further refined in 
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To start the MEM analysis of a given meson correlation function we need to specify a default 
model that incorporates all prior knowledge about the spectral function we want to determine. 
Before presenting the analysis based on default models related to the fits discussed in the previous 
subsections we performed a MEM analysis in analogy to what we had done earlier [3| ■ We used as 
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FIG. 9. The ratio of second and zeroth thermal moment of the correlation function GairT) obtained from 
fits with different values for the continuum cut-off parameter uio/T and fixed A^^. Circles show results 
of a MEM analysis where the fits have been used as default model. The band gives the result extracted 
from the continuum extrapolated correlation function. The lower curve shows the corresponding free field 
(infinite temperature) value, which is about 5% smaller than the value obtained at T" ~ 1.45Tc. 




FIG. 10. Spectral functions obtained from a maximum entropy analysis using the free spectral function as 
default model. 



default model the free spectral function, i. e. an ansatz that docs not include a contribution linear in 
uj at small energies but rather pfl^^{u!)/uj ~ w^/T. This default model together with the resulting 
output spectral function is shown in Fig. (TU] Apparently this ansatz also forces p""^'"(cj)/a; to 
vanish at w = 0. In particular, it leads to a suppression of the dilepton rate at small energies. The 
small contribution at small energies is compensated by large contributions at a; ~ (2 — A)T. This 
feature has been observed also in our first analysis of vector spectral functions 18 1. 

To allow a MEM analysis to find a contribution p{u}) ~ w at the origin this possibility needs to 
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FIG. 11. Spectral functions obtained from a maximum entropy analysis using for the default model the 
spectral functions shown in Fig. [Sfleft). The right hand figure shows the difference between the output 
spectral function obtained from the MEM analysis and the input spectral function used in each case. 

be offered by the default modej^ It thus seems to be appropriate to use the result of the analysis 
presented in the previous subsections as default model in a MEM analysis. It definitely summarizes 
the best knowledge we have at present about the spectral function in the vector channel. Using the 
un-modified version of the Breit-Wigner plus continuum ansatz, Eg. IV. 21 as default model changes 
results only little and leads to an output spectral function that is compatible with the input form. 
It also reproduces the result obtained for the electrical conductivity given in Eq. IV. 91 

In order to judge the stability of the spectral function obtained with our fit ansatz, we used a 
class of default models of the form given in Eq. IV. 101 We performed a MEM analysis using as 
input default models the spectral functions shown in Fig. |8l We controlled the error on the output 
spectral functions by performing the MEM analysis on jackknife blocks. In this way we find that 
the error of the output spectral functions never is larger than 10% for small energies and is below 
5% for uj/T>8. In Fig. [TTJlcft) we show results for spectral functions obtained in a MEM analysis 
by using a set of default models with varying luq and fixed A^^. The right hand figure shows the 
difference between the input default models and the output spectral functions. We note that these 
differences correspond to changes in the spectral functions which are generally smaller than 5%. 
They increase in particular for small values of w as the x^/d.o./. of the input default model gets 
worse. We also note that the MEM analysis reproduces the calculated ratio of thermal moments 
as well or even better than the fits used as a default model, although the value of the thermal 
moments itself did not directly enter the MEM analysis. 

VI. ELECTRICAL CONDUCTIVITY, THERMAL DILEPTON AND PHOTON RATES 

Here we want to discuss the consequences of our results for thermal dilepton rates as well as the 
electrical conductivity and the related soft photon rate. 



On the other hand, in correlation functions for quantum number channels where no transport contribution is 
present, MEM is flexible enough to suppress or eliminate completely a contribution p(A) ~ A included in the 
default model. This is, for instance, the case in the pseudo-scalar channel. 
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FIG. 12. Thermal dilepton rate in 2-flavor QCD (left). Shown are results from fits without a cut-off on the 
continuum contribution (loq/T = 0) and with the largest cut-off tolerable in our fit ansatz {ljq/T = 1.5). 
The HTL curve is for a thermal quark mass itit/T = 1 and the Born rate is obtained by using the 
free spectral function. The right hand part of the figure shows the spectral functions that entered the 
calculation of the dilepton rate. 



In Fig. [12] we show the thermal dilepton rate calculated from Eg. III.14l for two massless {u, d) 
flavors. We use the results obtained with our Breit-Wigner plus continuum fit ansatz, Eq. IV. 21 
as well as results obtained with a truncated continuum term. For the latter we use the case, 
ljq/T = 1.5, Aui/T = 0.5, which gave a x^/d-o.f of about 1. These results are compared to 



a dilepton spectrum calculated within the hard thermal loop approximation 12| using a thermal 
quark mass ttit/T = 1. Obviously the results are in good agreement for all U!/T>2. For l<cj/r<2 
differences between the HTL spectral function and our numerical results is about a factor two, 
which also is the intrinsic uncertainty in our spectral analysis. At energies ui/T<\ the HTL results 
grow too rapidly, as is well known. 

In the limit w — > the results for ph{lo)/uj, and thus also for the electrical conductivity, are 
sensitive to the choice of fit ansatz. Within the class of ansatzc used by us a small value of 
Pii{uj)/ui seems to be favored. Our current analysis suggests, 

2 < lim ^^^^ < 6 at T ~ 1.45 . (VI.l) 

This translates into an estimate for the electrical conductivity 

1/3 < < 1 at T~ 1.45 . (VL2) 

Using Eq. IIL15l this yields for the zero energy limit of the thermal photon rat^, 
H R 

lim oj—r^ = (0.0004 - 0.0013) Tr* ~ (1 - 3) • 10"'^ GeV^ at T ~ 1.45 . (VI.3) 



Here we used Tc ~ 165 MeV. This is a value relevant for QCD v^ith 2 light quarks rather than the critical 
temperature for a pure SU (3) gauge theory. 
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VII. CONCLUSIONS 

At a fixed value of the temperature, T ~ 1.45Tc, we liave performed a detailed analysis of vector 
correlation functions in the high temperature phase of quenched QCD. A systematic analysis at 
different values of the lattice cut-off combined with an analysis of finite volume and quark mass 
effects allowed us to extract the vector correlation function in the continuum limit for a large 
interval of Euclidean times, 0.2 < tT < 0.5. In this interval the correlation function has been 
determined to better than 1% accuracy. Furthermore, we determined its curvature at the mid- 
point of the finite temperature Euclidean time interval, tT ~ 1/2. 

We analyzed the continuum extrapolated vector correlation functions using several fit ansatze 
that differ in their low momentum structure. We find that the vector correlation function is best 
fitted by a simple ansatz that is proportional to a free spectral function plus a Breit-Wigner term. 
Already this simple three parameter ansatz yields a good description of the correlator with a 
small x^/d-o.f.. Other fits with a /d.o.f. < 1 are possible and suggest that the low energy 
structure of the spectral function currently has a systematic uncertainty of about a factor two. 
Some generic features of the spectral function are, however, robust. For energies U!/T>{2 — 4) the 
spectral function is close to the free form. At present we used an ansatz proportional to the free 
spectral function to parametrize the spectral function in this regime. However, the hard thermal 
loop spectral function should give a good description at these energies and can be incorporated in 
our ansatz, which cuts off the small energy continuum contributions. 

For small energies, uj/T<{1 — 2), the spectral function is significantly larger than the free result, 
but smaller than the HTL spectral function, which diverges at small energies; at energies uj/T ~ 1 
the thermal dilepton rate is about an order of magnitude larger than the leading order Born 
rate. In order to analyze quantitatively to what extent this enhancement can account for the 
experimentally observed enhancement of dilepton rates at low energies 0, \^ we will need results 
on the spectral function at temperatures closer to the transition temperature as well as knowledge 
on its momentum dependence. With this a complete analysis of dilepton rates that takes into 
account the hydrodynamic expansion of dense matter created in heavy ion collision will become 
possible j43|. 
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